The central question is: how should we allocate a fixed budget to three PrEP support interventions to maximize the number of HIV infections averted over the course of 10 years among MSM in Atlanta?
The model parameters representing policy levers being considered:
• POIP – percentage of MSM reached by initiation intervention
• POAC – adherence intervention capacity (max. number of people who can receive the intervention per year)
• PORC – retention intervention capacity (max. number of people who can be active in the intervention in any given week)
These are synonymous with “Uptake Capacity”, “Adherence Capacity”, and “Retention Capacity”.
The outcome we are looking at in Figure 1 is the average proportion of active PrEP users who are highly adherent. We expect this proportion to go up with increasing adherence intervention capacity. At a certain point, this highly adherent proportion stops increasing with adherence capacity because all new PrEP users are already enrolled in the adherence intervention. When the uptake intervention is well-funded (High PrEP Uptake), a greater adherence capacity can be utilized. Consequently, we see the proportion of active PrEP users who are highly adherent start to plateau at a lower adherence capacity when PrEP uptake is also low. “Low PrEP Uptake” and “High PrEP Uptake” are coded as the simulations with the lowest and highest quartiles of the parameter that controls the initiation intervention capacity.
The y-axis on Figure 2 denotes the mean proportion of the PrEP eligible population that is actively taking PrEP for each simulation over the course of the 10 year time horizon. Increasing the rate of PrEP uptake and the average length of time each PrEP user stays on the medication will both raise overall coverage levels. Increasing retention intervention capacity raises coverage levels until the capacity exceeds the number of PrEP users who can utilize the capacity. Increasing PrEP uptake shifts this point of excessive retention capacity to a higher value. Visually speaking, it appears that marginal effect of retention capacity on PrEP coverage is the same at all levels of PrEP uptake, except in the space where retention capacity exceeds demand. If infections averted is a function of PrEP coverage, then this suggests that the lack of a synergistic interaction between the uptake and retention interventions when retention demand is not exceeded. “Low PrEP Uptake” and “High PrEP Uptake” are coded as the simulations with the lowest and highest quartiles of the parameter that controls the initiation intervention capacity.
Figure 3 shows the same dynamics in a slightly different way than Figure 2. The parameter represented by the x-axis on Figure 3 (PrEP Intervention Initiation Percentage) is the percentage of individuals newly indicated for PrEP who become actively engaged in the initiation intervention. The higher retention capacity values do not lead to an increase in PrEP coverage, and the relationship between coverage levels and the initiation intervention capacity parameter appears to be roughly linear and independent of retention capacity.
So far, we have only described the relationship between the parameters controlling the capacity of PrEP support interventions and intermediate outcomes in an exploratory fashion. To answer our research question within a mathematical optimization framework, we need to reduce the the messiness of the simulation data down to smooth functions of the policy parameters of interest.
We use a generalized additive model (GAM) to fit a flexible and smooth function to our simulation output. In this statistical modeling approach, our outcome is the percent of cumulative infections averted over the course of 10 years, and the predictors of interest are uptake program initiation proportion, retention program capacity, and adherence program capacity. Like generalized linear models, GAMs are specified by a link function and a mean variance relationship (family). Because our outcome (percent infections averted) is positive and continuous, we used a Gamma family mean-variance relationship. Even though the percentage of cumulative infections averted is technically bounded between 0 and 100, we never see enough infections averted for this upper bound impact our choice of family for the statistical model. We chose to use a log link function because it provided the best fit to our data in terms of the minimized generalized cross-validation (GCV) score, which is used for smoothness selection in the mgcv package in R. In GAMs, relationships between the individual predictors and the dependent variable follow smooth patterns that can be nonlinear, which is necessary for capturing the relationships that we observed in the exploratory analysis of the simulation output.
We can visualize the generalized additive model predictions for infections averted to get a sense of how the program parameters interact and affect outcomes. On the three dimensional and contour plots, the axes span the entire allowable range of each capacity parameter (Uptake Capacity: 0-0.10, Adherence Capacity: 0-2000, Retention Capacity: 0-2000).
The y-axis denotes the predicted number of infections averted based on the generalized additive model and three different intervention capacity values. The plot has three facets for the minimum, median, and maximum values for the uptake intervention parameter, and each curve is color-coded according the minimum, median, and maximum adherence intervention capacity. This allows one to visualize the effect of changing retention capacity on infections averted at various combinations of the other program capacities.
These plots are faceted by uptake intervention capacity and color coded according to retention intervention capacity. We are able to see how adjusting retention intervention capacity affects predicted infections averted at representative combinations of the other two intervention capacity parameters.
These plots are faceted by adherence intervention capacity, and color coded by retention intervention capacity. We are able to see how adjusting uptake intervention capacity affects predicted infections averted at representative combinations of the other two intervention capacity parameters.
The generalized additive model informs the leverage of each intervention capacity parameter in maximizing infections averted, but we also need to define the monetary costs of increasing the capacities for each intervention in order to translate the optimization problem into economic terms. Using micro-costing data and the resources described in the clinical trial protocols, we estimated constant costs of increasing program capacities by a single unit, which is defined differently for the three interventions. The resource costs and personnel costs that go into operating the interventions are higher than the costs used here, because only certain costs would be need to be covered by the CDC budget. The perspective of this analysis is that of the CDC, and therefore we have excluded all costs that would ultimately be paid for by sources such as insurance.
The parameter that reflects the capacity of the uptake intervention, POIP, controls the percentage of MSM who are newly indicated by PrEP that are reached by uptake intervention. Once reached by the intervention that individual’s probability of starting PrEP is close to 100%. In order for 1% of all newly indicated MSM to be reached by the intervention within the simulation, we make the assumption that 1% of the entire potential population for the intervention (undiagnosed MSM who are not already on PrEP) must be already using the app. In the EpiModel simulation of around 12,000 individuals, roughly 8,000 are undiagnosed MSM not on PrEP. A POIP value of 1% implies that 80 of these individuals are actively using the app, and we calculate advertising costs based on what it would take to reach that many individuals. In estimating advertising costs, we make the following assumptions: -500 ad clicks are required to convert a single active app user -$0.79 per click (based on Facebook advertising costs) -App users will use the app for an average of 6 months -Advertising costs must be sustained over the full 10 years of the time-horizon
Additionally, PrEP eligibility will be determined by the app based upon HIV tests and STI tests, among other things. While the tests themselves will not have to be paid for through the CDC budget***, it is reasonable to expect that the shipment and handling of these tests must be covered by the CDC. It is assumed that the CDC will have to pay for the shipment of 2 HIV testing kits ($7.00 each) and 2 STI testing kits ($16.25 each) per year for each individual who is actively using the app.
***Now assuming 10% cost-sharing for HIV/STI tests
It is expected that the behavioral counseling sessions associated with the Life-Steps intervention will not be covered by insurance as they are not normally required for PrEP use, and therefore the bulk of the costs must be covered by the CDC. Each individual who starts the intervention is assumed to engage in 6 1-hour long counseling sessions. The personnel costs for a registered nurse paid at the median hourly wage in Atlanta for 6 total hours of work is $315.18. The overhead of utilizing clinic space was calculated using the CMS outpatient provider fee schedule and comes to $246.36. The parameter in the simulation controls the number of individuals who are allowed to start the adherence intervention each year, and because the intervention is run for 10 years, the unit cost associated with increasing this yearly capacity by 1 slot must be multiplied by 10.
The parameter controlling retention capacity works a little differently than the parameter controlling adherence capacity in the simulation. Whereas the adherence capacity parameter defines the number of individuals who are allowed to start the adherence intervention in a year, the retention capacity parameter (PORC) defines the number of individuals who are allowed to be actively utilizing PrEP@Home at any given time. Therefore, each unit of retention capacity represents an increase of 10 person-years of utilization of the services provided by PrEP@Home. While on the intervention the shipment of 4 HIV testing kits ($7.00 each) and 2 STI testing kits ($16.25 each) must be paid for by the CDC per person-year of intervention utilization. Additionally, it is assumed that, on average, 15 minutes of call-in support will be required per person-year which translates to a cost of $12.44 per person-year.
The other services associated with PrEP@Home such as the diagnostic tests themselves and telehealth appointments are far more costly, but because they are required in standard PrEP care, it is expected that health insurance rather than the CDC will have to pay ***
***Now assuming 10% cost-sharing for HIV/STI tests
In summary, the following linear unit costs are assumed for the budget optimization problem, based on the assumptions outlined above:
$69,130,560.00 per percentage point increase to PrEP Optimization Initiation Proportion (POIP)
$5,650.80 per yearly adherence intervention capacity slot (POAC)
$1,104.64 per retention intervention capacity slot (PORC)
With the assumed unit costs and the model-derived predictions for infections averted as a function of the intervention capacities, we use a nonlinear optimization algorithm to solve for the optimal allocation of funds at a range of budgets from $700,000 to $6,000,000.
In Figure 13, we are looking only at how our predicted percentage of 10-year cumulative infections averted increases as a function of the budget constraint. At 100 different total budget constraints (shown on the x-axis), the optimization algorithm yields a set of capacity parameter values and the associated outcome predicted by the GAM model. Without seeing how the budget is allocated, Figure 13 tells us how much overall predicted value is extracted from a given budget. At the lower budget values, the marginal increase in infections averted is greater than at higher budget values.
In Figure 14, we visualize the overall cost per infection across the range of budget values being considered. It is apparent that higher budget values exhibit slightly diminished returns in terms of their cost-efficiency in averting infections.
Rather than assessing the overall cost per infection averted over the range of budget values, Figure 15 considers the average marginal increase cost per infection averted across the range of budget values. This can be conceptualized as the rate at which cost efficiency is worsening with diminished returns and other interactions.
Figures 16, 17, and 18 show the optimal capacity parameter values at a range of budget constraints. Note that the range of allowable adherence capacity is 0-2000, and the noise around an adherence capacity of 0 in Figure 16 means that the adherence intervention never represented a cost-efficient use of funds at any budget value. At the end of this document, we will look a sensitivity analysis in which we find the cutoff point in the unit cost of the adherence intervention at which the intervention becomes viable to some extent in the budget range being considered, because it clearly does not not enter into the realm of possibilities under the current assumptions.
Figure 17 tells us that at low budget values, it is not optimal to allocate funding to the uptake intervention. Once the retention intervention starts to give significantly diminished returns, injecting money into the uptake intervention becomes optimal. Once the uptake/initiation intervention starts to become funded around a $1M budget, the level of funding appears to increase roughly linearly as the budget increases.
Figure 18 shows that the retention intervention is the most cost-efficient and absorbs the entire budget at low budget values. Around a retention capacity of 1,350, the population of active PrEP users has been effectively saturated in its engagement with the retention intervention. The only way to extract significant additional value from the retention intervention is to also increase the number of individuals who are on PrEP, and the uptake intervention provides a means for achieving this. As the uptake intervention raises the size of the PrEP active population, it opens the door for additional cost-efficient funding of the retention intervention.
Figure 19 displays the fraction of the budget that is allocated to each intervention capacity as a function of the total budget size. At low budgets, the retention intervention is the most utilized lever for averting cumulative infections, but as the allowed budget increases, spending on the uptake intervention begins to dwarf spending on the retention intervention.
Figure 20 shows the same general relationship as Figure 19, except the y-axis is the total amount of funding rather than the fraction of funding allocated to each intervention capacity across various budget sizes. As the uptake intervention becomes more prominent, it does not displace any of the money previously allocated to the retention intervention. The most cost-efficient use of the budget is always to provide retention intervention capacity to individuals on PrEP who can utilize it. However, once the entire population of active PrEP users is utilizing the retention intervention, there is no additional utility associated with increasing its capacity. The funding for the retention intervention rises slightly at higher budgets because the initiation intervention causes more people to start PrEP, which raises the ceiling for usable retention capacity.
In Table 1 we can see the numerical summaries for the economic and epidemiological outcomes associated with a selection of the budgets considered. One takeaway is that the cost per infection averted is over three times higher at the highest budget value ($6 million) compared to the lowest budget value ($700,000).
Table 1:
The adherence intervention does not receive any investment in the optimized budget allocations derived from our estimated costs per unit capacity for the three interventions. Through a sensitivity analysis, we found that the adherence intervention receives significant funding across a range of budget constraints when the cost per unit capacity is 48% of our estimated value. Figure 21 shows the optimal fraction of the total budget that is consumed by each intervention as a function of total budget size under the new cost assumption for the adherence intervention. Even when the cost per unit capacity for the adherence intervention is reduced by 52%, the retention intervention still absorbs the entirety of the funding at the lowest budget values. Once the retention intervention has been funded extensively enough to exhibit significantly diminishing returns, the adherence intervention is funded instead of the uptake intervention. Once the adherence intervention exhibits diminishing returns as well, further funding is allocated to the uptake intervention.
Figure 22 shows the optimal funding amount (in total $) allocated to each intervention as a function of total budget size.